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INTRODUCTION 


The primary goal of the Laser Goodynamics Satellite 
Mission is to employ precision laser tracking observations 
to measure the dynamic behavior of the earth at the 
centimeter level of accuracyt The currently available 
precision of the laser observations is betto'' than five 
centimeters and is progressively improving. LAGEOS is in a 
more stable and predictable orbit than other retro- 
reflector-carrying satellites, and the perturbation model 
for the satellite's motion is improving with each sot of new 
observations. Possible errors in the laser observations 
themselves will become increasingly Important in their 
application to derive important geodetic parameters such as 
the relative movement and deformation of the earth's 
tectonic plates. 

Many of the errors encountered in a laser ranging 
system can be reduced only through the employment of 
Improved instrumentation. Thsy arise primarily in the time 
interval measurement resolution and precision, the ability 
to keep real time, and the stability of the system 
oscillator. If wo exclude from consideration satellite- 
dependent error sources such ns the size or depth of the 
target and the effect of coherent fading of the laser beam 
duo to retro-reflector non-planarity, the most significant 
remaining source of error is the correction of the ranging 
measurement for atmospheric refraction. 
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The range ref met ion correction la auaceptlblo to 
errors in our assumptions concerning atmospheric compoaltlon 
and homogeneity, ns well as in the numerical approximation 
to an analytical model. A further source of error in the 
refraction model arises in the measurements of the 
meteorological conditions on which the approximations to the 
atmosphere are based. These quantities are limited by 
instrumentation aecurncy and can never be completely free of 
operator observation and transcription errors unless they 
are automatically recorded. 

We have considered the refraction effect from three 
perspectives. An analysis of the axioms on which the 
currently accepted correction algorithms wore b«, sod vjas the 
first priority. The integrity of the meteorological 
measurements on which the correction model is based was also 
considered and a large quantity of laser observations was 
processed in an effort to detect any serious anomalies in 
them. The effect of refraction errors on geodetic 
parameters estimated from laser data using the most recent 
analysis procedures was the focus of the third element of 
our stnvlies. The reported resnltvS concentrate on refraction 
errors which vjc have found to he critical in the eventual 
use of the data for measurements of crustal dynamics. 
Petaila of analyses in which refraction or data error was 
found to be insignificant are not reported. One of the 
criteria which we have adopted to determine significance is 
that a satellite range observation taken at 20® elevation or 
a tower calibration measurement of 10 km. two-way distance 
be perturbed by more than one centimeter* Several of our 
experlmentvS merely confirmed the results published by other 



workers in the field of atmospheric refraction, 
occasionally dwelt upon the details of individual 
measurements to Illustrate points concerning the 
conclusions which we discuss in the final section 
report. 


We have 
sets of 
general 
of this 
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REFRACTIVITY AT OPTICAL FREQUENCIES 


The determination of the range between a laser system 
and the target presumes a correct value of the velocity of 
propagation of the laser energy as well as the time it takes 
the laser energy to travel to the target and return. In 
propagating through the atmosphere, the laser energy is 
slowed and suffers geometric bending. Both effects can be 
explicitly determined if the true index of refraction along 
the path of propagation is known. For ranging systems where 
accuracy and precision to a few tens of centimeters is 
required, simple atmospheric models and nominal measurements 
of atmospheric conditions are sufficient. However, for 
those systems where range accuracy better than a few 
centimeters is desired, more complex atmospheric models and 
better atmospheric parameters are required for accurate 
correction of the propagation velocity. 

The atmospheric refraction correction applied to GSFC 
laser range data is computed from the formulation developed 
by Marini and Murray (Ref. 0). The correction is computed 
as a function of the laser wavelength, the station 
coordinates, the local temperature, pressure and relative 
humidity, and the elevation angle of the satellite. Safety 
requirements restrict the elevation angle at which the 
lasers acquire observations to be above twenty degrees. 

The Dry Term 

The formula adopted by Marini and Murray for the phase 
refractlvity N of dry air is that established in Resolution 
No. 1 of the 13th General Assembly of the l.A.G. (Ref. 1): 
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ORIGINAL PAGE 13 
OF POOR QUALITY 


N « 10^ (n-1) 


287.604 + 

X2 x"* 


* I * 273.15 
T 1013.25 


. .43908 . .00367 

77.53174 + T + 7 — 


* 


P 

T 


in which n is the refractive index of dry air 

X is the wavelength of light 

P is the atmospheric pressure in mbars 

and T is the temperature in °K. 

The expression is based upon the work of Barrell and 
Sears (Ref. 2) who investigated the refraction and 
dispersion of dry, C 02 -free air by means of an interference 
ref ractometer . Th^^y give an accuracy of jH). 01*10”^ to their 
formula and quote a term (also adopted in the liarini and 
Murray treatment) for the effect of humidity. 

The group refractivity Ng of air is arrived at 
through application of the dispersion formula 
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m 

Quality 


Ng . N . 


X 

dX 


N - X 


-2*. 43908 4^,00367 


77.53174 + 4 . .J.Q. 

X^ X^ 


- 80.343 


0.9650 4 . ^228 

X2 


T 


- 80.343 f(X) 


« 80.343 for red light of X - .6943u 


and 


« 80.343*1.023 * •= for green light of X * .5322u 


Green light will thus be refracted by an amount 2.3^ 
greater than red light for given atmospheric conditions. 

In the comparison by liacrell and Sears of their 
formula with those given by earlier workers they note good 
agreement (within .04X) with Perard (Ref. 3) and Rosters and 
Lampe (Ref. 4). However, their results differ from those of 
Meggers and Peters (Ref. 5) by 0.21Z for red light and 0.27% 
for green light. Meggers and Peters would thus give a 
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refraction correction approximately 2 cm. shorter than the 
Marini Murray formulas for a 20° elevation range observation 
using a green laser for which the total refraction effect is 
about 7 meters. 

The important effect of Meggers and Peters' 
apparently erroneous dispersion formula on standard 
spectroscopic tables which had been in use for many years, 
was pointed out by Edlen (Ref. 6). In a later paper, Edlen 
(Ref. 7) compared more recent experimental results (up to 
1961) with his own dispersion formula for which he claimed 
an accuracy of 1 part in 10^ and which agrees with that due 
to Barrel! aud Sears to better than 1.4*10“® (.0005%). We 
must therefore conclude that the expression for group 
refractivlty chosen by Marini and Murray is accurate enough 
to describe the behavior of light in the visible spectrum to 
sub-millimeter levels for any practical laser ranging 
measurement . 

The Wet Terra 

The effect of water vapor in the atmosphere is given 
in the Marini Murray formulation as 


= -11.3 


_e 

T 


; in which 

w 

N is the wet term or g 'oup refractivlty 


3 
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and 


OF POOR QUALIT/ 


e is the partial pressure of water vapor in mbar 
which can 'j'je obtained from 


e 


Ry. 

* 6 . 11 * 10 ** 


7.5(T-273.15) 

237.3+(T-273.15) 

- 


v^here is the relative humidity. 

At a temperature of 300°K, a relative humidity of 
100% would yield a partial pressure e ■= 35 mbar which gives 
a wet component of refractivity amounting to 1.33*10"’^ (or 
about 0.5%). This could amount to 1.5 cm in the 3 meter 
round trip tower calibration correction for a tower located 
5 km from the laser site. To restrict the error due to 
refraction in the wet atmosphere to one centimeter in such a 
calibration measurement, the relative humidity only needs to 
be measured very crudely and need only be applied if it is 
greater than 50%. To limit the error in the dry 
refractivity term to one centimeter in the tower calibration 
observation, the average pressure over the tower distance 
must be known to 4 mbar and the average temperature to 1°C. 
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REFRACTION CORRECTION OF SATELLITE MEASUREMENTS 


The model adopted by Marini and Murray for the 
correction of an optical path through the vertical 
atmosphere was based on an exhaustive study of the 
literature on atmospheric correccion formulae for tracking 
data. To match the centimeter level accuracy of the laser 
systems at low elevation anglesu the integral evaluations of 
the group i ndcx of refraction along the phase path given by 
Saastamolnen (Ref, 8) are incorporated into Marini's 
continued fraction form of the range correction (Ref, 9). 

Marini and Murray found that the relative accuracy of 
the finally fchosen refraction model was better than one 
centimeter when compared to corrections computed from ray- 
tracing radiosonde profiles. They point out that errors 
caused by factors common to both methods are not in evidence 
in their tests. These include the equations for the group 
refractive index, the errors in which we have found to be 
small. The common assumption of hydrostatic equilibrium is 
noted by the authors, and the hydrostatic equation on which 
they base their atmospheric model is also implicit in the 
ray-tracing method through its use to infer the heights of 
the radiosonde observations. We are unable to question this 
axiom of atmospheric modelling and have in fact relied on 
the hydrostatic equation in a similar treatment of 
radiosonde observations which we describe below in an 
attempt to place a bound on a final source of error in the 
Marini Murray model: the assumption of horizontal 
homogeneity. 
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OF POOK “ 

refraction correction ZiR in meters to a laser 
range measurement finally arrived at in the liarini Murray 
trr. .traent is 


CR 


f (X) 
f (<j),H) 


si nE+ 


A + B 

B / ( A+B )~ 
8 inE+0 . 0 1 


in which 

A » .002357P + .0001416 

2 

B - ( 1.084*10~8) + (4.734* 10"®) * '(3-i ' /kT 

K » 1.163 “ 0.00968 cos 2<f) - 0.00104T + 0.00001435P 

E is the true elevation of the satellite 
f((|),H) = 1 - .0026 cos 2c|> - .00031H 

<j) is the latitude of the laser I'^ite 
H is the altitude of the laser site in km. 
and f(X)j P, T, and e are as defined above 
The expression can be linearized to the approximation 
AR ~ .0072P + 0.0004e + .00000025PT 

for E “ 20 ^ 
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and 


P » 


lOOOffibar, T •* 300°K and e *■ 20mbar(50XR.H» ) 


If the error in the range correction is to hs limited 
“o one centimeter the surface pressure must therefore be 
known to l.babar, the water vapor pressure to within 60 
units of percentage relative humidl^y and the surface 
temperature to 40°C, These accuracies are clearly within 
the range of properly calibrated instruments and errors in 
the refraction correction due to problems with surface 
measurements will be due to any undetected variation in them 
over the satellite pass and errors in reading and recording 
the observations. 


Variations in the Real Atmosphere 

The real atmosphere varies both temporally and in 
three-dimensional space. The time-variation of atmospheric 
turbulence has bean found by Gardner (Ref. 10) to affect 
satellite range measurements to an insignificant level for 
most comb itia tions of turbulence strength, scale size and 
propagation path length. However, for horizontal paths near 
the earth's surface, such as those used for laser system 
delay calibration links, centimeter level errors were found 
due to turbulence. 

The effects of an atmosphere which varies in three 
dimensions have been investigated by Iyer and Bufton (Ref. 
11), who show that assymetry in the refractive index could 
be. accounted for by expressing the range correction as a 
series of terms. An estimate of the higher order correction 
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terme was obtained and can be seen to be In quantitative 
agreement with the horizontal gradient effect determined by 
Gardner (Ref. 12). Meteorological data from the Haven Hop 
network of weather stations located near Washington, D.C. 
was considered in each of these studies. A correction 
formula which compensates for the gradient effect wa? 
developed by Gardner (Ref. 12) and evaluated by Gardner et 
al. (Ref. 13) by ray tracing through J)-D refractivity 
profiles generated using the Project Haven Hop radiosonde 
measurements . 

Gardner et al. conclude that the observed errors in 
the refractivity observations can be reduced using the 
gradient correction formula of Ref. 12. They note that, at 
a 20° elevation angle the sea level gradient error is 
approximately 5mm. for a horizontal surface temperature 
gradient of l°C/100kra. The increase in refractivity caused 
by the negative temperature gradient from the thermal 
equator toward the colder climates at the poles was clearly 
evident in the data and gave rise to the dominant systematic 
component of the gradient correction formula, in which the 
most significant effect on a range in meters is given by 

X + C 08 C y) V(PTK) 


in which 


^(X),E,P,T and K are as defined above, 
V is the gradient operator, 
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X is the east unit vector, 
y Is the north unit vector, 
and o is the satellite azimuth angle. 

A laser range observation taken at 20® elevation 
passes through the atmosphere at lOkra, altitude at a 
distance of 10 cot 20® “ 27 km. from the laser site. If the 
effects on a range measurement of variations in atmospheric 
conditions up to 10 km. altitude are to be considered, 
observations at the spacing of a few tens of kilometers 
should be studied. In the following section we describe an 
assessment of the effect on the refraction of laser ranges 
of horizontal gradients measured at this short spatial 
interval . 


Analysis of Radiosonde Observations in New Mexico 

A series of radiosonde observations was collected 
during the Prototype Artillery Subsystem (PASS) Project 
(Ref. 14). One of the purposes of this project was to 
collect a large volume of atmospheric data to be used for 
further research and development, particularly in the areas 
of sound ranging applications and ballistics. Observations 
of temperature, pressure and relative humidity were 
collected in November 1974 by radiosonde instruments mounted 
on balloons launched from several sites in White Sands, New 
Mexico. The close temporal separation of many of the 
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vertical profile measurements (within one half hour) 
provided a useful attribute for studies of horir.ontal 
gradient effects. The close spatial separation of the 
launch sites (several dozen kilometers) allows an assessment 
of shorter wavelength variations than were observable in the 
Haven Hop Project, in which the stations were separated by 
several hundred kilometers. 

On the other hand, the limited vertical profile 
extent (up to 300 mbars) placed a restriction on the 
measured vertical profiles, and the close spacing of the 
sites limited the grid-size required for detecting the 
trends at lower spatial frequencies that were provided by 
the Haven Hop observations. 


Three Dimensional Ray-Tracing Procedure 

A three-dimensional ray tracing code was developed to 
investigate real and theoretical variations in horizontal 
ref ractivity , The three dimensional r-afractive structure of 
the atmosphere is specified by interpolation among vertical 
refractivity profiles specified above (up to) nine grid 
points on the ground. The ray path linking the laser and 
the satellite is obtained by choosing an approximate initial 
ray direction and integrating the ray trace equation 


ds 


/ dr. 


Vn 
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outwards toward the satellite. Here, n, denotes the phase 
refractive index, s is a distance along the ray path, r is a 
vector designating a point on the ray, and 7 is the gradient 
operator. The direction and amount of miss is used to 
update or correct the initial ray direction and a new path 
is computed. This procedure is iterated until the ray hits 
the satellite center with an acceptably small error. 

Once the correct ray path has been found, the 
refraction correction is computed and consists of two 
terras: the group delay correction 


/(n '’l)ds 
8 


and a geometric correction 


/ ds-R. 


In these equations, n denotes the group refractive index 

O 

and R is the true satellite range along a straight line 
path. 
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Numoricnl Methods 


The ref ractivlty n(h)-l and its vertical 
gradient 9i\/3h are constructed above iiach of the grid points 
on the ground at regular intervals of height. The 
refractivity is computed using Owens formulae (Ref. IS) from 
vertical profiles of temperature, pressure and relative 
humidity. When the refractivity at an arbitrary altitude h 
and horizontal position (X,Y) is required, it is obtained by 
interpolation. First, values of u(h) and 9n/Dh are obtained 
by four point hagr«nge interpolation on each of the vertical 
profiles. N(U,X,Y), 9n( h , X , Y) / 8h , 3n/9G and 9n/9(j) 
(where 0 and <|) are latitude and longitude) are then obtained 
from a linear two dimensional interpolation using up to six 
grid points. 

The integration of the ray path differential equation 
is curried out using a modified Hamming prod Ic tor‘-corrector 
method with variable stop size, initiated by an iterated 
Runge-Kutta scheme. Once the ray path has been found (using 
a Regain Falsi technique for iteration), the group delay is 
computed. The path is divided Into segments which are short 
near the ground and increase in length using a twelve point 
Gaussian quadrature and summed to give the total group delay 
correction. 

The temperature profile assumed for altitudes above 
the radiosonde data was that given in the U.S. Standard 
Atmosphere 1976 (Ref. 16) for altitudes at ll, 20 , 32 , 47 , 
51, 71 and 84.5 km. This standard profile is shown in 
Figure I and was used to supplement the radiosonde data 
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above 300mbar (about 9 km.) and alternatively to extrapolate 
upwards a profile based only on surface temperature and 
assuming a constant lapse rate up to the 11 km. value of the 
Standard Atmosphere profile. 

The surface observation of pressure was then combined 
with the temperature profile through the hydrostatic 
equation to produce the corresponding pressure profile. 


Assessment of the Vertical Profile Observations 

The radiosonde observations were of limited value for 
making a full assessment of spatial trends in the atmosphere 
due to the limited spacing of the launch sites and the 
drifting of the balloons during ascent. The coordinates of 
balloons launched within a one half hour interval on 
November 2nd, 197A is shown in Figure 2. By the time the 
balloons had reached 500 mbar (about 7 km.) they had drifted 
through a distance comparable to the launch site spacing. 
However it can be seen from Figure 2 that they had 
maintained approximately the same horizontal separation 
during ascent, and can therefore be used to monitor 
horizontal gradients in a consistent space whose upward 
direction is at a nearly constant angle to the vertical. 

In order to simplify the quantitative assessment of 
the radiosonde observations to which we are limited, the 
launch sites were placed on a regular rectangular grid 
(denoted as "ideal position" in Figure 2) with a longitude 
spacing of 24 km. end a latitude spacing of 16 km. Our 


17 


model assumed a laser site L placed at the grid point shown, 
tracking at a number of elevation angles through the 
profiles above the grid points and at azimuths in the 
directions N, NW, W, SW and S of L. The ideal positions of 
the launch sites are always farther from the laser than 
their real positions, and thus the horizontal gradient 
levels inferred in the ray- tracing procedure will be 
conjervative. 

The real temperature profiles for the experiment of 
November 2nd are shown in Figure 3, together with those 
approximated by a constant lapse rate between the surface 
measurements and the 11 km. value of the Standard Atmosphere 
profile. A common feature between 4 and 5 km. in each 
profile Indicates the capability of the observations to 
monitor real variations. The shaded areas of Figure 3 are 
measures of temperature changes which would not be 
adequately reflected in the constant lapse rate profiles. 

The differences in temperature from the laser profile 
as a function of altitude for the November 2nd data is shown 
in Figure 4, which indicates a temperature range of about 
5°C at each altitude level. Although the observation error 
cannot be excluded as a source of these apparent temperature 
gradients, the overall consistency of the trends suggest 
that variations of a few degrees are present in the real 
atmosphere. The precision of radiosonde instrument has been 
assessed at 1°C at the surface varying linearly to 2.5°C at 
30 km. (Ref. 17). Temperature differences of several 
degrees at the surface and at altitude are shown in Figure 5 
for the November 2nd experiment together with two other 
cases a few days later. 
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The Effect of Temperature Gradients on Refractivity 

Under the assumptions of vertical ascent and linear 
extrapolation outside the ideal grid, the ray-tracing 
procedure predicts the effects on a range measurement at 20° 
elevation as shown in Figure 6a. Figure 6b indicates the 
refraction effect if a constant temperature lapse rate from 
the surface to 11 km. is assumed. The level of the effects 
is similar for either model of temperature profile but there 
is little similarity in. the patterns predicted by each 
model . 


In Figure 7 the refraction effect on a 20° elevation 
range observation on November 2nd is shown for the models 
described above, together with that due to Marini-Murray , 
which assumes no horizontal gradient between the profiles. 
The contribution to the refraction effect from each segment 
of the atmospheric profiles is shown in Figure 7 and plotted 
in Figure 8 as the difference from the Marini-Murray model 
value. The effect of temperature gradients on refraction is 
seen to be greatest at about 10 km. The very large 
refraction effect at this altitude predicted by the real 
data balances a contribution of opposite sign at lower 
altitudes. 

The detailed analysis of real temperature 
observations indicates very high sensitivity of the 
refraction effect to the temperature gradients inferred by 
our procedure and is complicated by the drifting of the 
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balloons. The assumption of a constant lapse rate below 11 
km. based on surface temperature yields a more stable 
refractivity effect and corresponds more closely to the 
assumptions of Marini-Murray and Gardner (Ref. 12) about the 
behavior of the atmosphere. We therefore adopted this 
simpler model in a more exhaustive study of the full PASS 
data set. 

All cases in which there were surface observations of 
temperature and pressure within one half hour from each of 
the six sites shown in Figure 2 were chosen. The ideal site 
positions were adopted and the effect on a 20® elevation 
range from the laser site L was computed for azimuths 
corresponding to the N, NW, W, SW and S directions. The 
effect in centimeters is shown in Figure 9 for each of the 
thirty- tliree available cases. The table of Figure 9 is 
divided into four time periods of approximately one week 
each and the average weekly effect tabulated for each 
direction in Figure 10. Neither the weekly breakdown of the 
refraction effect nor the total variation over the 
approximately monthly time span suggests systematic 
variation in azimuth over the 24 km by 16 kra spacing. The 
variation can be seen to amount to almost 2 cm. over the 
full experimental period. The implications of these results 
are discussed in the las'fc section of this report. 
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INTEORiTY «F METKOROLOGIGAL OBSERVATIONS 


During tlie early raontha of 1979 four mobile laaera 
wore co-located with the GSFG stationary system (STALAS) in 
Greenbolt, Maryland. The relative locations of the aystoaa 
are shown at the top of Figure lla and they are denoted by 
the letters A, 8, G and 1). The heights of the laaera on a 
reference ellipsoid, determined from the laser observations 
taken during the co-location testa were STALAS? i5 m., A? 
'5 m., B; 14 m. , G? 14 ra. and D? 6 m. The meteorological 
measurements collected by each system to drive the 
refraction model were available in the format for the 
ranging observations and were tabulated for each pass of 
LAGEOS data during which at least one mobile laser obtained 
observations simultaneously with STALAS. 

The time of day, the value of the temperature, 
pressure and relative humidity at STALAS and the difference 
between the simultaneously obsorviug stations' 
raoteorologlcnl reading and that of STALAS are tabulated in 
Figure lla for three occasions in February 1979. Figures 
11b and lie show similar information for simultaneous LAGEOS 
pass acquisitions in March and April 1979 respectively. The 
most striking difference in the readings at each site is in 
relative humidity which fortunately is the parameter to 
which the refraction correction is least sensitive. 
However, errors in the relative humidity as large as 59S (on 
3“22 at 10 hrs) would produce errors of almost one centimeer 
in the refraction correction of a 20'^^ elevation range or in 
that of a calibration measurement at 10 km two-way tower 
distance. Ther-^j appears to bo no systematic difference in 
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relative humidity, pressure or temperature between the 
various sites but the difference between the sites A, B and 
STALAS, which are located within a few dozen meters can be 
seen to amount to and 4 rabars on occasion. It must be 

assumed that this difference is due to errors in the reading 
of the Instruments for such closely spaced observations. An 
error of 4^C would produce an error of 4 cm in a 10 km two- 
way tower calibration measurement and would certainly 
confound any attempt to monitor horizontal temperature 
gradients to apply a gradient correction term. A pressure 
error of 4 mbar would cause a refraction error of 0.4% in 
all satellite ranges taken from a station at sea level. 
This would amount to a little less than 3 cm for a 20^ 
elevation ranging observation. 

In Figures 12a, b and c the differences in 
temperature and pressure radlngs taken at the stationary 
Cireenbelt laser (STA) and a mobile laser (MOB) are shown, 
(or a collocation period in early 1980. Although 

temperature differences as large as 3°G are occasionally 
recorded, no systematic temperature difference can be 
seen.% On the other hand, the pressure recorded at the 
mobile site is systematically lower than Chat at the 
stationary laser. The average bias of the mobile system 
during the months of March and April amounts to -3.5 mbar. 
In a comparison of these measurements with those collected 
by the stations of the National Weather Service, Gibbs and 
Mayer (Ref. 18) suggest that the error occurred in the 
barometer at the mobile site. Large pressure biases were 
also found at the Haystack, Mass, and Patrick Air Force 
Base, Florida stations as well as large, random temperature 
differences with the NWS stations in Goldstone, California. 



GEODETIC PARAMETER ESTIMATION 


The approach conventionally adopted for the 
estimation of geodetic parameters from laser range data is 
based on a Bayesian weighted least squares orbit 
determination scheme such as the GEODYN system (Ref. 19). 
Laser data, preprocessed to Include calibration and 
atmospheric refraction corrections, are reduced to 
simultaneously estimate the satellite orbit, components of 
station location at one or more sites, and possibly orbit 
model characteristics, such as an atmospheric drag or solar 
radiation pressure scaling factor. 

The accuracy of geodetic quantities such as relative 
station location estimated from laser observations is 
critically dependent on the design of the experiment. When 
data from near-earth satellites are used, the analyst's main 
preoccupation is to design an experiment to reduce the 
effects of dynamic force model error. The final data 
configuration has usually been that in which as many 
observations have been acquired in as short a period of time 
as possible. More subtle designs will be required to reduce 
the effect of refraction model error. 

A simple test has been made, based on a technique 
used to measure relative station heights and interstation 
chord distances between lasers in the West Atlantic tracking 
GEOS-3 (Ref. 20). Single pastes of data of less than ten 
minutes in length were employed to reduce the effect of the 
dominant source of error; the geopotential model. The 
geodetic measurements were made with a precision of 
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approximately 15 cm. if all known error sources arc 

considered and a perfect refraction model is assumed. If we 
assume the rather large error of 1% in the refractive 
modulus used in preprocessing the observations, the inter- 
station baselines Increase by up to 8 centimeters. The 
orbital fits to the range data remained at the observation 
noise level of 7 cm. even when the refractive modulus was 
Increased by 10%, as the refraction error was completely 
absorbed by the estimated parameters. Refraction errors 
will actually arise in a more random pattern than the 
elevation-dependent bias assumed in this test, which 
nevertheless suggests that the refraction model must be 
carefully considered as an error source in geodetic 

parameter estimation at the centimeter level. 

Error Analysis of TLRS Tracking LAGEOS 

The LAGEOS orbit is high enough that it is much less 
sensitive to drag, solar radiation pressure and geopotential 
model error than are other laser geodetic satellites in 
near-earth orbits. The accuracy of currently available 
dynamic orbit models allows orbital fits to the range data 
at the noise level of less than 10 cm. for several 
revolutions and of less than 25 c-m. for orbital arc lengths 
of several days (Ref. 21). At this level of precision the 
effects of instrumentation and refraction error may approach 
that in the dynamic model. 

To investigate the effect of refraction error in 
comparison to other sources of error, a LAGEOS tracking 
configuration in the western United States suggested by 
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Christodoulidis and Smith (Reft 22) was considered# It is 
depicted graphically in Figure 13 and is comprised of two 
stations considered fixed at San Diego (SANDIE) and at 
Quincy, California# The relative positions of transportable 
lasers at T1 and T2 (a north-south baseline) or at T1 and T3 
(an east-west baseline) were the parameters whose estimates 
were investigated using the ORAN error analysis system (Ref. 
23) under a variety of circumstances# 

For the case in which two transportable lasars were 
available for simultaneous occupation of the sites, all 
possible LAGEOS ranging measurements above 20° elevation 
angle from the four stations were simulated during a 5 day 
time period. The noise level of the observations was 
assumed to be 10 cm. at a repetition rate of 1 second# The 
error analysis system computes the effects of perturbations 
to the force and measurement model on the estimated 
parameters, which in this case were the six elements of the 
orbit and three position components of each transportable 
system. For the cases in which 100% efficiency was assumed 
for the laser systems, the effects on the baseline between 
the transportable systems caused by errors in the listed 
unadjusted parameters are shown in Figure 14 in the columns 
labelled 100%. The two alternative configurations (north- 
south and east-west) are represented graphically above the 
corresponding table. The error sources considered were 25 
cm„ in each component of the location of the fixed stations, 
10 cm. biases in the laser range observations, a 1% error in 
ref ractivity , an error in the universal constant GM of 1 
part in 10^ and a measure of gravity error in the variance- 
covariance matrix of the GEM-9 geopotential model (Ref. 
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24). The total error in the north-south baseline is shown 
in Figure 14 to be 2.9 cm. and contains a significant effect 
due to refraction and instrumentation biases. These 
observation model errors are less important in the east-west 
baseline error estimate, which amounts to 2.2 cm., largely 
due to GM error. A more realistic situation in which 50% 
efficiency is assumed for the laser systems was also 
simulated and the results shown in Figure 14o In this case 
only 50% of all available passes were considered in the data 
reduction scheme, A high sensitivity to geopotential model 
error is indicated in the baseline error estimates which 
amount to 11.3 cm. for the north-south configuration and 6.2 
cm. east-west. Tne effects of refraction and 
instrumentation bias do not significantly Increase when the 
system efficiency is decreased to 50%, 

A situation in which we are afforded the luxury of 
only a single transportable laser ranging system was also 
considered. In this case the transportable system must 
occupy the sites between which the baselines must be 
estimated in separate time periods. In the simulations, 
consecutive 5-day occupations were assumed and Figure 15 
shows the results for 100% and 50% values of system 
efficiency. The degradation in system efficiency is seen to 
be far more critical in the case of consecutive 
transportable laser site occupations than in the case where 
simultaneous ranging is possible. In particular the effect 
of 1% refraction error amounts to 3.9 cm. in the east-west 
baseline error estimate. In order to approach the results 
possible with 5 days of simultaneous data, the lasers 
require 30 days of continuous tracking at 50% efficiency in 
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two coneGcutlve site occupations of a single transportable 
systotn. The effects of the unraodelled errors on baseline 
estimates from this extended period of deployment is also 
shown in Figure 15. The effects of refraction error are 

reduced to 1,6 cm, in the north-south baseline and 0.7 cm. 
in the east-west bnsellno estimate. 
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DISCUSSION 


The Marlni~Murray refraction model which Is applied 
to the laser data collected as part of the Crustal Dynamics 
Program Is appropriate within the limitations stated by its 
authors. Marini and Murray's treatment is based on some 
simplifying assumptions about the meteorological conditions 
above the laser site. It was found to produce results in 
essential agreement with those given by our ray tracing 
procedure, even when constraints on the variations in 
vertical profiles of temperviture were relaxed. The 
correction algorithm is simple and the coi'rection is well- 
defined from easily observed surface measurements. 

The main limitation of the current model is in its 
requirement for spherical symmetry in the atmosphere above 
the laser site. This restriction is stated by the authors 
and has been studied by several workers since the original 
correction algorithm was published. The observations upon 
which studies of horizontal gradients had been based before 
our analysis were made at meteorological observing stations 
whose closest separation was about 100 km. and whose largest 
separation was about 600 km. (see Ref. 13). This data 
configuration revealed relatively small average effects on 
the refraction of laser range observations. The stations 
were spaced far enough apart to detect the refraction effect 
due to the temperature gradient between the pole and the 
equator. 


A laser observation of 20° elevation passes through 
the lower tropopause at a horizontal distance of 30 or 40 
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kilometers from the station. It is at this altitude that 
any horizontal gradients in the atmosphere would be most 
critical. We have considered possible gradients as measured 
by radiosonde balloons separated by distances between 10 km, 
and 35 km. The benefits of a spacing commensurate with the 
surface distance travelled by a low elevation ranging signal 
are however limited by the errors in the extrapolation 
necessary to infer the atmospheric structure outside our 6- 
station grid. The short time intervals (less than half an 
hour) between observations of the vertical profile were 
compatible with the time span of a satellite pass and would 
therefore indicate the effects of any gradients such a pass 
might experience. The limited altitude (300 mbars) of the 
vertical profile measurements prevented any direct 
measurement of temperature gradients closer to the 
tropopause. Our results therefore suffer from possible 
errors of extrapolation in the vertical direction, which 
would tend to exaggerate any real variation. 

Some stability was restored to the gradient 
measurements by assuming a model for conditions aloft 
similar to that on which the Mar ini-Murray model was 
based: a constant temperature lapse rate based on surface 
observations. Errors in horizontal extrapolation still 
remain with this approach and will exaggerate gradients 
measured in a small grid, in which there is little 
redundancy of information to eliminate effects with very 
short spaclal wavelength or those due to instrumental 
error. The simplification of the vertical atmospheric 
structure also eliminates the effect of balloon drift which 
made the detection of any systematic trends in the gradients 


29 


difficult from the real vertical profile observations. A 
range measurement at 20° elevation passing through an 
atmosphere with the computed horizontal gradients would 
diffur by several centimeters from that in a spherically 
stratified atmosphere. 

The analysis of the real observations was finally 
reduced to the consideration of qualitative evidence of 
temperature gradients at altitude. Systematic patterns of 
temperature variations were detected in vertical profiles 
which were evenly spaced in the vertical and relative 
horizontal directions, as they were measured from balloons 
with nearly consistent drift rates. These temperature 
variations (up to 5°C at surface and at altitude) are too 
large to be included in the accepted error budget for the 
instruments. Their systematic nature supports our 

contention that real gradients may exist and are not 
artifacts of measurement error. The inference we draw that 
horizontal gradients may contribute to several centimeters 
of error in a pass of low elevation ranging observations is 
not incompatible with the results of others who have 
averaged over larger space and time intervals and found 

smaller effects. We feel however that we have erred in the 
direction of exaggeration and that these results be 

considered upper bounds on the effects of refractlvity 
variation. 

Assessments of the integrity of the surface 
observations which are collected at a laser station for 

inclusion in the refraction correction were not 

reassuring. Variations of several degrees in temperature 
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were found between observations made at collocated laser 
sites for which identical meteorological conditions would be 
expected. This variation could be the result of operator 
error in reading and recording the instruments or occasional 
instrument malfunction. It could also be caused by natural 
variation of temperature over a distance of a few hundred 
meters. Holdahl (Ref. 25) has observed temperature 

variations of several degrees centigrade due to 
stratification in the boundary layer close to the ground. 
Temperature error in the currently adopted refraction 

correction procedure will mostly affect tower calibration 
measurements which are made by the currently deployed laser 
systems to allow for cable delays which would bias the 

satellite measurements. Internal calibration procedures 
would eliminate this possible error source, which could 
however be reduced if more temperature observations were 

made in the vicinity of the laser stations and particularly 
in the direction of the calibration tower. Regular 

calibration of each station's barometer would reduce the 
possibility of pressure errors, which dominate the 
refraction error budget for satellite range corrections. 

The precision of the pressure and temperature 
readings reported with the laser ranging observations 
slightly complicated our attempt to confirm that the Marini~ 
Murray model had in fact been applied to the ranges. The 
reported values are truncated to the nearest millibar and 
degree, although the refraction correction is computed using 
more precise observations. However, in light of the 
physical and possibly instrumental limitations on the 
meteorological values and their application, we do not 
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consider this a serious source of error. Of all the 
perturbations to the refraction model which we have 
investiRated , this is the simplest to eliminate. 

Tl\e evidence of variability in meteorological 
conditions which wo have found is not completely 
unexpected. Variations in pressure and temperature are 
extremely difficult to predict or model. Wo feel that the 
elimination of consequent errors in the refraction model is 
more effectively accomplished by designing geodetic 
parameter estimation procedures which include redundancy. 
The simplest way to introduce thlvS element is to extend the 
period over which satellite observations are collected to 
yield a single geodetic measurement. Simulation studios 
based on a typical experiment have Indicated that larg© 
refraction errors Cvsn he considerably reduced by extending 
the campaign duration from *5 days to 30 days. This time 
period is currently required to also eliminate errors in 
geodetic measurements caused by lack of knowledge of the 
satellite perturbation model and the relative location of 
supporting stations. Those other elements of tlve error 
budget will improve with time duo to expected improvements 
in instrumentation and in the development of a better 
satellite force model. Improvements in the atmospheric 
model are unlikely to keep pace with progress in other areas 
of satellite laser analysis and we must therefore bo careful 
to maintain the. Integrity of the meteorological observations 
on which our refraction model depends. 
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FIGURE 2. BALLOON COORDINATES AT 880 MBAR AND 500 MBAR (11-2-10^1 
(SHOWS DRIFT TO NORTH-EAST) 
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FIGURE 3 REAL TEMPERATURE PROFILES AND PROFILES FROM SURFACE DATA 

(THE DIFFERENCE IS SHADED) 
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figure 6b. EFFECT ON 20° ELEVATION RANGE DUE TO GRADIENTS COMPUTED FROM SURFACE DATA 
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FIGURE 11b. METEOROLOGICAL OBSERVATIONS DURING 1979 COLLOCATION (MARCH) 
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FIGURE 12a. METEOROLOGICAL CONDITIONS DURING 1980 COLLOCATION (FEB) 
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FIGURE 12b. METEOROLOGICAL CONDITIONS DURING 1980 COLLOCATION (MARCH) 
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FIGURE 12c. METEOROLOGICAL CONDITIONS DURING 1980 COLLOCATION (APRIL) 
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